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Abstract 

An exact re-normalization group equation describes the dependence of the free energy 
on an infrared cutoff for the quantum or thermal fluctuations. It interpolates between 
the microphysical laws and the complex macroscopic phenomena. We present a simple 
unified description of critical phenomena for 0(iV)-symmetric scalar models in two, three 
or four dimensions, including essential scaling for the Kosterlitz-Thouless transition. 
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1 Introduction 

Particle physics and statistical physics share a common theoretical challenge. For given 
"microphysical laws" observable "macrophysical" quantities at much larger characteristic 
length scales have to be computed. This concerns, for example, the computation of the hadron 
spectrum from perturbative QCD or the thermodynamic equation state for a liquid of helium 
from the atomic interactions. The difficulties arise from the importance of fluctuations. More 
technically, the "transition to macrophysics" requires the solution of a complicated functional 
integral. Beyond the identical conceptual setting of particle and statistical physics there is 
also quantitative agreement for certain questions. The critical exponents at the endpoint of 
the critical line of the electroweak phase transition (the onset of crossover for the critical mass 
of the Higgs scalar) are believed to be precisely the same as for the liquid-gas transition near 
the critical pressure. This reflects universality. 

Exact renormalization group equations aim for a piecewise solution of the fluctuation 
problem. They describe the scale dependence of some type of "effective action". In this 
context an effective action is a functional of fields from which the physical properties at a 
given length or momentum scale can be computed. The exact equations can be derived as 
formal identities from the functional integral which defines the theory. They are cast in the 
form of functional differential equations. Different versions of exact renormalization group 
equations have already a long history The investigation of the generic features of their 

solutions has led to a deep insight into the nature of renormalizability. In particle physics, the 
discussion of perturbative solutions has led to a new proof of perturbative renormalizability 
of the 4 -theory || . Nevertheless, the application of exact renormalization group methods to 
non-perturbative situations has been hindered for a long time by the complexity of functional 
differential equations. Some type of expansion is needed if one wants to exploit the exactness 
of the functional differential equation in practice - otherwise any reasonable guess of a realistic 
renormalization group equation does as well. 

Since exact solutions to functional differential equations seem only possible for 
some limiting cases, it is crucial to find a formulation which permits non-perturbative 
approximations. These proceed by a truncation of the most general form of the effective 
action and therefore need a qualitative understanding of its properties. The formulation of 
an exact renormalization group equation based on the effective average action || [| has been 
proven successful in this respect. It is the basis of the non-perturbative flow equations which 
we discuss in this review. The effective average action is the generating functional of one- 
particle irreducible correlation functions in the presence of an infrared cutoff scale k. Only 
fluctuations with momenta larger than k are included in its computation. For A; — ^ all 
fluctuations are included and one obtains the usual effective action from which appropriate 
masses and vertices can be read off directly. The k dependence is described by an exact 
renormalization group equation that closely resembles a renormalization group improved one- 
loop equation |§. In fact, the transition from classical propagators and vertices to effective 
propagators and vertices transforms the one- loop expression into an exact result. This close 
connection to perturbation theory for which we have intuitive understanding is an important 
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key for devising non-perturbative approximations. Furthermore, the one-loop expression is 
manifestly infrared and ultraviolet finite and can be used directly in arbitrary number of 
dimensions, even in the presence of massless modes. 

The aim of this review is to show that this version of flow equations can be used in 
practice for the transition from microphysics to macrophysics, even in the presence of strong 
couplings and a large correlation length. It is largely based on a more extensive report by J. 
Berges, N. Tetradis, and the author |I0| . (We refer to this report as well as to ]TT[ for a more 
extensive list of references.) We derive the exact renormalization group equation and various 
non-perturbative truncations in section ^|. In section 3 we discuss the solutions to the flow 
equation in more detail for 0(A)-symmetric scalar models. We derive the scaling form of the 
flow equation for the effective potential and compute critical exponents for second-order phase 
transitions in three dimensions. Finally we discuss essential scaling for the Kosterlitz-Thouless 
phase transition for two-dimensional models with a continuous Abelian symmetry. 

Several results in statistical and particle physics have been obtained first with the method 
presented here. This includes the universal critical equation of state for spontaneous breaking 
of a continuous symmetry in Heisenberg models [ 12 1 , the universal critical equation of state for 
first-order phase transitions in matrix models []13| |, the non-universal critical amplitudes with 
an explicit connection of the critical behavior to microphysics (CO2) [Q, a quantitatively 
reliable estimate of the rate of spontaneous nucleation |15|, jlBf , a classification of all possible 
fixed points for (one component) scalar theories in two and three dimensions in case of a 
weak momentum dependence of the interactions ||17|| ) the second-order phase transition in 
the high temperature 4 quantum field theory [18|], the phase diagram for the Abelian Higgs 



model for N charged scalar fields |19[ gy]. In particle physics it has led to the prediction that 
the electroweak interactions become strong at high temperature, with the suggestion that 
the standard model may show a crossover instead of a phase transition |21|, [22|| . In strong 
interaction physics the interpolation between the low temperature chiral perturbation theory 
results and the high temperature critical behavior for two light quarks has been established 



in an effective model PB[. All these results are in the non-perturbative domain. In addition, 



the approach has been tested successfully through a comparison with known high precision 
results on critical exponents and universal critical amplitude ratios and effective couplings. 

Our main conclusion can be drawn already at this place: the method works in practice. 
If needed and with sufficient effort precise results can be obtained in many non-perturbative 
situations. New phenomena become accessible to a reliable analytical treatment. 



2 Non- Pert ur bat ive flow equation 

2.1 Average action 

We will concentrate on a flow equation which describes the scale dependence of the effective 
average action T^. ||. The latter is based on the quantum field theoretical concept of the 
effective action T, i.e. the generating functional of the Euclidean one-particle irreducible {IP I) 
correlation functions or proper vertices. For the vacuum or ground state this functional is 



3 



obtained after "integrating out" the quantum fluctuations. The masses of the excitations, 
scattering amplitudes and cross sections follow directly from an analytic continuation of the 
1PI correlation functions in a standard way Furthermore, the field equations derived from 
the effective action are exact as all quantum effects are included. In thermal and chemical 
equilibrium T includes in addition the thermal fluctuations and depends on the temperature T 
and chemical potential fi. In statistical physics T is related to the free energy as a functional 
of some space-dependent order parameter <p{x). For vanishing external fields the equilibrium 
state is given by the minimum of V. More generally, in the presence of (spatially varying) 
external fields or sources the equilibrium state obeys 

M J{x), (2.1) 



5<j>{x) 

and the precise relation to the thermodynamic potentials like the free energy F reads 

F = ™- + 0-tJ*°M*W*). (2.2) 

Here 4> eq (x) solves fl2.1|) , T eq = T[(f) eq ], and N is the conserved quantity to which the chemical 
potential is associated. For homogeneous J = j/T the equilibrium value of the order 
parameter is often also homogeneous. In this case the energy density e, entropy density s, 
"particle density" n and pressure p can be simply expressed in terms of the effective potential 
U{<t>) = TT/V, namely 

dU 8U dU j<f> 

€ = U - T dT-^> S = ~dT + T' 

n = -^,p = -U= -TT/V. (2.3) 

Here U has to be evaluated for the solution of dU/d<ft = j , n = N /V and V is the total volume 
of (three-dimensional) space. Evaluating U for arbitrary yields the equation of state in the 
presence of homogeneous magnetic fields or other appropriate sources. 

More formally, the effective action T follows from a Legendre transform of the logarithm of 
the partition function in the presence of external sources or fields (see below). Knowledge of F 
is in a sense equivalent to the "solution" of a theory. Therefore T is the macroscopic quantity 
on which we will concentrate. In particular, the effective potential U contains already a large 
part of the macroscopic information relevant for homogeneous states. We emphasize that 
the concept of the effective potential is valid universally for classical statistics and quantum 
statistics, or quantum field theory in thermal equilibrium or the vacuum^. 

The average action is a simple generalization of the effective action, with the distinction 
that only fluctuations with momenta q 2 > k 2 are included. This is achieved by implementing 
an infrared (IR) cutoff ~ k in the functional integral that defines the effective action T. In 

1 The only difference concerns the evaluation of the partition function ZorW = lnZ = —(F — /j,N)/T. For 
classical statistics it involves a D-dimensional functional integral, whereas for quantum statistics the dimension 
in the Matsubara formalism is D + 1. The vacuum in quantum field theory corresponds to T —> 0, with V/T 
the volume of Euclidean "spacetime" . 
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the language of statistical physics, Tk is a type of coarse-grained free energy with a coarse 
graining length scale ~ k~ l . As long as k remains large enough, the possible complicated 
effects of coherent long-distance fluctuations play no role and T k is close to the microscopic 
action. Lowering k results in a successive inclusion of fluctuations with momenta q 2 > k 2 and 
therefore permits to explore the theory on larger and larger length scales. The average action 
Tfc can be viewed as the effective action for averages of fields over a volume with size k~ d || 
and is similar in spirit to the action for block-spins on the sites of a coarse lattice. 

By definition, the average action equals the standard effective action for k = 0, i.e. T = T, 
as the IR cutoff is absent in this limit and all fluctuations are included. On the other hand, in 
a model with a physical ultraviolet (UV) cutoff A we can associate Ta with the microscopic 
or classical action S. No fluctuations with momenta below A are effectively included if the IR 
cutoff equals the UV cutoff. Thus the average action Tk has the important property that it 
interpolates between the classical action S and the effective action r as k is lowered from the 
ultraviolet cutoff A to zero: 

T A ^S , limr fc = r. (2.4) 

The ability to follow the evolution to k — > is equivalent to the ability to solve the theory. 
Most importantly, the dependence of the average action on the scale k is described by an exact 
non-perturbative flow equation presented in the next subsection. 

Let us consider the construction of I\ for a simple model with real scalar fields Xa, 
a = 1 ... N, in d Euclidean dimensions with classical action S. We start with the path 
integral representation of the generating functional for the connected correlation functions in 
the presence of an IR cutoff. It is given by the logarithm of the (grand) canonical partition 
function in the presence of inhomogeneous external fields or sources J a 

W k [J] = \nZ[J] =\nj D X exp {-S[ X } - AS k [ X ] + j d d xJ a (x) X a (x)^ . (2.5) 

In classical statistical physics S is related to the Hamiltonean H by S = H/T, so that e~ s 
is the usual Boltzmann factor. The functional integration j D\ stands for the sum over all 
microscopic states. In turn, the field Xa( x ) can represent a large variety of physical objects like 
a (mass-) density field (N = 1), a local magnetisation (N = 3) or a charged order parameter 
(N = 2). The only modification compared to the construction of the standard effective action 
is the addition of an IR cutoff term A<S&[x] which is quadratic in the fields and reads in 
momentum space (Xa(~ Q) = Xa(l)) 

ASk[x] = \j -^- d Rk ^a(-q)x a (q)- (2.6) 

Here the IR cutoff function R k is required to vanish for k — > and to diverge for k — > oo (or 
k — > A) and fixed q 2 . This can be achieved, for example, by the exponential form 

Rk(q) = If - (2.7) 
e q 1 - 1 
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where the wave function renormalization constant Z k will be defined later. For fluctuations 
with small momenta q 2 <C k 2 this cutoff behaves as R k (q) ~ k 2 and allows for a simple 
interpretation: Since AS/Jx] is quadratic in the fields, all Fourier modes of x with momenta 
smaller than k acquire an effective mass ~ k. This additional mass term acts as an effective 
IR cutoff for the low momentum modes. In contrast, for q 2 3> k 2 the function R k (q) vanishes 
so that the functional integration of the high momentum modes is not disturbed. The term 
AS k [x] added to the classical action is the main ingredient for the construction of an effective 
action that includes all fluctuations with momenta q 2 > k 2 , whereas fluctuations with q 2 < k 2 
are suppressed. 

The expectation value of Xi i- e - t ne macroscopic field 0, in the presence of ASfcfx] and J 
reads 

SW k [J] 



We note that the relation between and J is fc-dependent, 
In terms of W k the average action is defined via a modified Legendre transform 



(2.8) 



; [J] and therefore J = J k 



W k [J] + / d d xJ a (x)<f> a (x) - AS k [4>], 



(2.9) 



where we have subtracted the term ASk[<j)] in the rhs. This subtraction of the IR cutoff term as 
a function of the macroscopic field <fi is crucial for the definition of a reasonable coarse-grained 
free energy with the property r a ~ S. It guarantees that the only difference between T k 
and T is the effective infrared cutoff in the fluctuations. Furthermore, it has the consequence 
that T k does not need to be convex, whereas a pure Legendre transform is always convex by 
definition. The coarse-grained free energy has to become convex [23L only for k — > 0. These 
considerations are important for an understanding of spontaneous symmetry breaking and, in 
particular, for a discussion of nucleation in a first-order phase transition. 

In order to establish the property T\ m S we consider an integral equation for T k that 
is equivalent to (|2.9| ). In an obvious matrix notation, where Jx = / d d xJ a {x)x a {x) = 



I J0z J a(-p)Xa(p) and R k ,ab(<l,q') = R k (q)Sab(^) d S(q - q'), we represent (|2~5| ) as 



exv[W k [J) 



Dxexp 



As usual, we can invert the Legendre transform 



S[x] + Jx - -jXRkX 



to express 



J= 6 4^ + <pR k . 



(2.10) 



(2.11) 



It is now straightforward to insert the definition (|2.9|) into ( 2.10|) . After a variable substitution 
x' = X ~ one obtains the functional integral representation of 



exp(-r fc [0D 



Dx'exp [ -S[(f) + x] + -jtX 



\x' Rk x' 



(2.12) 



This expression resembles closely the background field formalism for the effective action which 
is modified only by the term ~ R k . For k — > oo the cutoff function R k diverges and the term 
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exp(— x'i? fc x'/2) behaves as a delta functional ~ 6\x'], thus leading to the property T k — > 5 
in this limit. For a model with a sharp UV cutoff A it is easy to enforce the identity T\ = S 
by choosing a cutoff function R k which diverges for k — > A, like -R^ ~ q 2 {e q2 l k2 — e q2 ^ 2 )~ l . 
We note, however, that the property Ta = 5 1 is not essential, as the short distance laws 
may be parameterized by T\ as well as by S. For momentum scales much smaller than A 
universality implies that the precise form of T\ is irrelevant, up to the values of a few relevant 
renormalized couplings. Furthermore, the microscopic action may be formulated on a lattice 
instead of continuous space and can involve even variables different from Xa{ x )- m this case one 
can still compute T\ in a first step by evaluating the functional integral (P-12|) approximately. 
Often a saddle point expansion will do, since no long-range fluctuations are involved in the 
transition from S to T\. In this report we will assume that the first step of the computation 
of Ta is done and consider Ta as the appropriate parametrization of the microscopic physical 
laws. Our aim is the computation of the effective action T from T\ - this step may be called 
"transition to complexity" and involves fluctuations on all scales. We emphasize that for large 
A the average action Ta can serve as a formulation of the microscopic laws also for situations 
where no physical cutoff is present, or where a momentum UV cutoff may even be in conflict 
with the symmetries, like the important case of gauge symmetries. 

A few properties of the effective average action are worth mentioning: 

1. All symmetries of the model which are respected by the IR cutoff AS k are automatically 
symmetries of r*.. In particular this concerns translation and rotation invariance. As a 
result the approach is not plagued by many of the problems encountered by a formulation 
of the block-spin action on a lattice. Nevertheless, our method is not restricted to 
continuous space. For a cubic lattice with lattice distance a the propagator only obeys 
the restricted lattice translation and rotation symmetries, e.g. a next neighbor interaction 
leads in momentum space to 

S= hl 7^E( 1 - cosa ^)x*(^) + --- (2-13) 
f 

The momentum cutoff < A, A = ir/a also reflects the lattice symmetry. A rotation 
and translation symmetric cutoff R k which only depends on q 2 obeys automatically all 
possible lattice symmetries. The only change compared to continuous space will be the 
reduced symmetry of IV 

2. In consequence, T k can be expanded in terms of invariants with respect to the symmetries 
with couplings depending on k. For the example of a scalar 0(N)-mode\ in continuous 
space one may use a derivative expansion (p = a a /2) 

T k = J d d x \u k (p) + l -Z^ k {p)d^ a d^ a + . . . j (2.14) 

and expand further in powers of p 

1_ 2 1_ 3 

Uk(p) = ^k{p-po{k)) + -7 fc (P - Po(k)) +■■■ 
Z^ k {p) = Zf tk (po) + Z'^ k {p ) (p - p ) + • • • (2.15) 
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Here p denotes the (fc-dependent) minimum of the effective average potential C4(p). 
We see that Tk describes infinitely many running couplings. 



3. Up to an overall scale factor, the limit k — > of Uk corresponds to the effective potential 
U = TT/V, from which the thermodynamic quantities can be derived for homogeneous 
situations according to eq. ( |2.3| ). The overall scale factor is fixed by dimensional 
considerations. Whereas the dimension of £4 is (mass) d the dimension of U in eq. 
(|2.3|) is (mass) 4 (for ft = c = ks = 1)- For classical statistics in d = 3 dimensions one 
has Uk = Tk/V and U = TUm^o Uk- For two-dimensional systems an additional factor 
~ mass appears, since Uk = T^/V^ = LTk/V implies U = TL~ l lirn&_>o Uk- Here L is the 
typical thickness of the two-dimensional layers in a physical system. In the following we 
will often omit these scale factors. 

4. The functional l\[0] = T %[(/>] + ASk[4>] is the Legendre transform of Wk and therefore 
convex. This implies that all eigenvalues of the matrix of second functional derivatives 
r^ 2 ) + Rk are positive semi-definite. In particular, one finds for a homogeneous field 4> a 
and q 2 = the simple exact bounds for all k and p 

U' k (p) > -R k (0) 
U'k(p) + 2pU' k \ P ) > -R k (0), (2.16) 

where primes denote derivatives with respect to p. Even though the potential U{4>) 
becomes convex for k — * it may exhibit a minimum at po(k) > for all k > 0. 
Spontaneous breaking of the 0(iV)-symmetry is characterized by limk^opo(k) > 0. 

5. The results for physical quantities have to be independent of the choice of the cutoff 
scheme Rk- On the other hand, both I\ and the flow with k are scheme- dependent. 
The scheme independence of the final results is a good check for approximations 



M, 27, 28, 29 . 



6. There is no problem incorporating chiral fermions, since a chirally invariant cutoff Rk can 
be formulated ^0, 31]. Gauge theories can be formulated along similar linesQ |32[ 33. 21 



BUB even though ASk may not be gauge invariant^. In this case the usual Ward 
identities receive corrections^ for which one can derive closed expressions [[3^, |36| . These 
corrections vanish for k — > 0. On the other hand they appear as "counterterms" in T\ 
and are crucial for preserving the gauge invariance of physical quantities. 

7. Despite a similar spirit and many analogies, there is a conceptual difference to the 
Wilsonian effective action . The Wilsonian effective action describes a set of different 
actions (parameterized by A) for one and the same model — the n-point functions are 
independent of A and have to be computed from by further functional integration. 



2 See also |4l] for applications to gravity. 

3 For a manifestly gauge invariant formulation in terms of Wilson loops see ref . [[58| . 

4 The background field identity derived in p2[ is equivalent p9|, Bol] to the modified Ward identity. 
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In contrast, Tk can be viewed as the effective action for a set of different "models" - 
for any value of k the effective average action is related to the generating functional of 
1PI n-point functions for a model with a different action Sk = S + AS^. The n-point 
functions depend on k. The Wilsonian effective action does not generate the 1PI Green 



functions 42 



2.2 Exact flow equation 

The dependence of the average action Tk on the coarse graining scale k is described by an 
exact non-perturbative flow equation |9|, |43|, f|5| 



1 d u 

dk Rk 



(2.17) 



The trace involves an integration over momenta or coordinates as well as a summation over 
internal indices. In momentum space it reads Tr = ^ a J d d q/(27r) d , as appropriate for the 
unit matrix 1 = (2n) d 5(q — q')S a b. The exact flow equation describes the scale dependence of 
Tk in terms of the inverse average propagator T k , given by the second functional derivative 
of Tk with respect to the field components 



,(2) 



ab 



S 2 T, 



50 a (-g)50 6 (g') 



(2.18) 



It has a simple graphical expression as a one-loop equation 



dTk 
dk 




with the full fc-dependent propagator associated to the propagator line and the dot denoting 
the insertion dkRk- 



Because of the appearance of the exact propag ator (T k 2) + R k )-\ eq. (plTp is a functional 



differential equation. It is remarkable that the transition from the classical propagator in 
the presence of the infrared cutoff, (S^ + Rk)^ 1 , to the full propagator turns the one-loop 
expression into an exact identity which incorporates effects of arbitrarily high loop order as 
well as genuinely non-perturbative effects^ like instantons in QCD. 

The exact flow equation ( |2.17| ) can be derived in a straightforward way ||. Let us write 



where, according to ( f2i) 



W fc [J] + / d d x.J(x)(t)(x) 



(2.19) 



(2.20) 



5 We note that anomalies which arise from topological obstructions in the functional measure manifest 
themselves already in the microscopic action T\. The long-distance non-perturbative effects ("large-size 



instantons") are, however, completely described by the flow equation (2.17). 
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and J = J k {<p). We consider for simplicity a one-component field and derive first the scale 
dependence of T: 



With <f)(x) = 5Wk/5J(x) the last two terms in (|2.21 ) cancel. The fc-derivative of W k is 



obtained from its defining functional integral ( |2.5| ). Since only R k depends on k this yields 

= (J^tx]) = (\J d d xd d y X (x)^R k (x,y) X (y)^ , (2.22) 
where R k (x,y) = R k (x — y) is the Fourier transform of R k (q) and 

(A[ X }) = Z- 1 J D X A[ X ] ex V (-S[ X } - A k S[ X ] + / d d xJ(x) X (x)). (2.23) 

Let G(x,y) = 5 2 W k /SJ(x)5J(y) denote the connected 2-point function and decompose 

(X(x)x(y)) = G(x, y) + ( X (x))(x(y)) = G(x, y) + (j>{x)<j>{y) . (2.24) 
After plugging this decomposition into (|2.22|) , the scale dependence of T k can be expressed as 



J^ t [0] = \J d d xd d y{^R k (x,y)G(y,x) + 4>{x)^R k (x,y)4>{y) 



The exact flow equation for the average action T k follows now through (|2.19 ) 



d „ r i„ r a „ l i 



SA; L ^ J 2 <% 2 



r fc [0] = ^Tr<{G--i4!> = ^Tr<' ri 2) [0] + i? fc ^ }► . (2.26) 



-i d 



dk' 



For the last equation we have used that T^\x,y) = 5 2 T k /5(J)(x)5(J)(y) = 5J(x)/5(j)(y) is the 
inverse of G(x,y) = 5 2 W k /5J(x)5J(y) = 5<p(x) / 5J(y): 

J d d yG(x, y) (ri 2) + R k ) (y, z) = S(x - z) (2.27) 

It is straightforward to write the above identities in momentum space and to generalize them 
to iV components by using the matrix notation introduced above. 
Let us point out a few properties of the exact flow equation: 

1. For a scaling form of the evolution equation and a formulation closer to the usual (3- 
functions, one may replace the partial /c-derivative in (|2.17|) by a partial derivative with 
respect to the logarithmic variable t = ln(/c/A). 
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2. Exact flow equations for n-point functions can be easily obtained from ( |2.17| ) by 

k 



(2) 

differentiation. The flow equation for the two-point function T k involves the three 



and four-point functions, 1^ and T^, respectively. One may write schematically 
d d 2 T k 



at k 



-i Tr 

2 



dR k d 



dt 



rjj. 2) + Rk 



-i 



r 



(3) 



rl 2) + R k 



Tr 



2 



dR k 



dt 

dRk 
dt 



rjjf + R k 



,(3) 



ri 2) + R k 



<3) 



-1 



rf + Rk 



r[ 2) + R k 



r 



(4) 



(2.28) 



By evaluating this equation for = 0, one sees immediately the contributions to the 
flow of the two-point function from diagrams with three- and four-point vertices. The 
diagramatics is closely linked to the perturbative graphs. In general, the flow equation 
for r[ involves and r[ n+2 - ) . 



3. As already mentioned, the flow equation (|2.17|) closely resembles a one-loop equation. 
Replacing 1^ by the second functional derivative of the classical action, , one obtains 
the corresponding one-loop result. Indeed, the one-loop formula for T k reads 



r fc [0] = S[0] + ^Trln (S^[<f>]+R k ) 



(2.29) 



and taking a fc-derivative of ( |2.29| ) gives a one-loop flow equation very similar to ( [2. 171 ). 



The "full renormalization group improvement" — > T k turns the one-loop flow 
equation into an exact non-perturbative flow equation. Replacing the propagator and 

(2) 

vertices appearing in T k by the ones derived from the classical action, but with running 
fc-dependent couplings, and expanding the result to lowest non-trivial order in the 
coupling constants, one recovers standard renormalization group improved one-loop 
perturbation theory. 



4. The additional cutoff function Rk with a form like the one given in eq. ( |2.15|) renders 
the momentum integration implied in the trace of ( |2.17| ) both infrared and ultraviolet 
finite. In particular, for q 2 k 2 one has an additional mass-like term Rk ~ k 2 in the 
inverse average propagator. This makes the formulation suitable for dealing with theories 
that are plagued with infrared problems in perturbation theory. For example, the flow 
equation can be used in three dimensions in the phase with spontaneous symmetry 
breaking despite the existence of massless Goldstone bosons for iV > 1. We recall that 
all eigenvalues of the matrix r^ 2 ^ + Rk must be positive semi-definite (cf. eq. ( |2.16|) ). We 
note that the derivation of the exact flow equation does not depend on the particular 
choice of the cutoff function. Ultraviolet finiteness, however, is related to a fast decay of 
dtRk for q 2 ^> k 2 . If for some other choice of Rk the rhs of the flow equation would not 
remain ultraviolet finite this would indicate that the high momentum modes have not 
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yet been integrated out completely in the computation of IV Unless stated otherwise 
we will always assume a sufficiently fast decaying choice of Rk in the following. 



5. Since no infinities appear in the flow equation, one may "forget" its origin in a functional 
integral. Indeed, for a given choice of the cutoff function all microscopic physics is 
encoded in the microscopic effective action IV The model is completely specified by 



the flow equation ( 2.17 ) and the "initial value" I\. In a quantum field theoretical sense 
the flow equation defines a regularization scheme. The "ERGE"-scheme is specified by 
the flow equation, the choice of Rk and the "initial condition" I\. This is particularly 
important for gauge theories where other regularizations in four dimensions and in the 
presence of chiral fermions are difficult to construct. For gauge theories I\ has to 
obey appropriately modified Ward identities. In the context of perturbation theory 
a first proposal for how to regularize gauge theories by use of flow equations can be 
found in We note that in contrast to previous versions of exact renormalization 

group equations there is no need in the present formulation to construct an ultraviolet 
momentum cutoff - a task known to be very difficult in non-Abelian gauge theories. 



6. Extensions of the flow equations to gauge fields |33], ^TJ, []34}|-[ji0| and fermions |30], |31 
are available. 



7. We emphasize that the flow equation ( [2.1 7|) is mathematically equivalent to the 



Wilsonian exact renormalization group equation 0, |4], ||, ||, [5J. The latter describes 
how the Wilsonian effective action changes with an ultraviolet cutoff A. Polchinski's 
continuum version of the Wilsonian flow equation can be transformed into eq. ( p. 17] ) 



by means of a Legendre transform, a suitable field redefinition and the association 



A = k j|3], |46|, |47| • Although the formal relation is simple, the practical calculation 
of Sjy from Tfc (and vice versa) can be quite involved^. In the presence of massless 
particles the Legendre transform of Tk does not remain local and is a comparatively 
complicated object. We will argue below that the crucial step for a practical use of 
the flow equation in a non-perturbative context is the ability to device a reasonable 
approximation scheme or truncation. It is in this context that the close resemblence of 
eq. ( p,17[ ) to a perturbative expression is of great value. 



In contrast to the Wilsonian effective action no information about the short-distance 
physics is effectively lost as k is lowered. Indeed, the effective average action for fields 
with high momenta q 2 ^> k 2 is already very close to the effective action. Therefore I\ 
generates quite accurately the vertices with high external momenta. More precisely, this 
is the case whenever the external momenta act effectively as an independent "physical" 
IR cutoff in the flow equation for the vertex. There is then only a minor difference 
between 1^ and the exact vertex T^ n \ 



6 If this problem could be solved, one would be able to construct an UV momentum cutoff which preserves 
^auge invariance by starting from the Ward identities for IV 
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An exact equation of the type ( |2.17| ) can be derived whenever R k multiplies a term 



quadratic in the fields, cf. ( |2.6| ). The feature that R k acts as a good infrared cutoff is 
not essential for this. In particular, one can easily write down an exact equation for the 
dependence of the effective action on the chemical potential flBfl . Another interesting 



exact equation describes the effect of a variation of the microscopic mass term for a 
field, as, for example, the current quark mass in QCD. In some cases an additional UV- 
regularization may be necessary since the UV-finiteness of the momentum integral in 
( |2.17 ) may not be automatic. 



2.3 Truncations 

Even though intuitively simple, the replacement of the (RG improved) classical propagator by 
the full propagator turns the solution of the flow equation (|2.17| ) into a difficult mathematical 
problem: The evolution equation is a functional differential equation. Once is expanded in 
terms of invariants (e.g. Eqs.( [2.14T ), ( |2.15|) ) this is equivalent to an infinite system of coupled 



non-linear partial differential equations. General methods for the solution of functional 
differential equations are not developed very much. They are restricted mainly to iterative 
procedures that can be applied once some small expansion parameter is identified. They cover 
usual perturbation theory in the case of a small coupling, the l/A^-expansion or expansions in 
the dimensionality 4 — d or d — 2. For the flow equation ( [2.17| ) they have also been extended 
to less familiar expansions like a derivative expansion |^9[ which is related in critical three- 
dimensional scalar theories to a small anomalous dimension ||50|| . In the absence of a clearly 



identified small parameter one needs to truncate the most general form of in order to reduce 
the infinite system of coupled differential equations to a (numerically) manageable size. This 
truncation is crucial. It is at this level that approximations have to be made and, as for all 
non-perturbative analytical methods, they are often not easy to control. 

The challenge for non-perturbative systems like critical phenomena in statistical physics or 
low momentum QCD is to find flow equations which (a) incorporate all the relevant dynamics 
so that neglected effects make only small changes, and (b) remain of manageable size. The 
difficulty with the first task is a reliable estimate of the error. For the second task the 
main limitation is a practical restriction for numerical solutions of differential equations to 
functions depending only on a small number of variables. The existence of an exact functional 
differential flow equation is a very useful starting point and guide for this task. At this point 
the precise form of the exact flow equation is quite important. Furthermore, it can be used 
for systematic expansions through enlargement of the truncation and for an error estimate in 
this way. Nevertheless, this is not all. Usually, physical insight into a model is necessary in 
order to device a useful non-perturbative truncation! 

Several approaches to non-perturbative truncations have been explored so far: 

(i) Derivative expansion. We can expand in the number of derivatives (p = ^(p a (f) a ) 

T k [</>] = J d d x \u k (p) + ^Z k (p)d^ a d^ a + \Y k (p)d,pd»p + 0(«9 4 ) j . (2.30) 
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The lowest level only includes the scalar potential and a standard kinetic term. The first 
correction includes the p-dependent wave function renormalizations Z^p) and Yk(p). 
The next level involves then invariants with four derivatives etc. 

One may wonder if a derivative expansion has any chance to account for the relevant 
physics of critical phenomena, in a situation where we know that the critical propagator 
is non-analytic in the momentum 7 . The reason why it can work is that the nonanalyticity 
builds up only gradually as k — > 0. At the critical temperature a typical qualitative form 
of the inverse average propagator is 

ri 2) ~gV + cfc 2 )-"/ 2 (2.31) 

with T) the anomalous dimension. Thus the behavior for q 2 — > is completely regular as 
long as k ^ 0. In addition, the contribution of fluctuations with small momenta q 2 <C k 2 
to the flow equation is suppressed by the IR cutoff Rk- For q 2 3> k 2 the "nonanalyticity" 
of the propagator is already manifest. The contribution of this region to the momentum 



integral in (|2.17|) is, however, strongly suppressed by the derivative dkRk- For cutoff 
functions of the type Q2.7D only a small momentum range centered around q 2 ~ k 2 
contributes substantially to the momentum integral in the flow equation. This suggests 
the use of a hybrid derivative expansion where the momentum dependence of r^— J d d xllk 
is expanded around q 2 = k 2 . Nevertheless, because of the qualitative behavior Q2.31| ), 
also an expansion around q 2 = should yield valid results. We will see in section |3| that 



the first order in the derivative expansion (|2.30|) gives a quite accurate description of 



critical phenomena in three-dimensional O(N) models, except for an (expected) error in 
the anomalous dimension. 

(ii) Expansion in powers of the fields. As an alternative ordering principle one may expand 
Tk in n-point functions 1^ 

rfc M=f;^ / (n^^o-^ijr^^,...,^). (2.32) 

n=0 ' \j=0 J 

If one chooses ||f] O as the fc-dependent expectation value of <fi, the series ( |2.32j ) starts 
effectively at n — 2. The flow equations for the 1PI Green functions 1^ are obtained 
by functional differentiation of ( |2.17| ). Similar equations have been discussed first in || 
from a somewhat different viewpoint. They can also be interpreted as a differential form 
of Schwinger-Dyson equations [^!| . 

(iii) Expansion in the canonical dimension. We can classify the couplings according to their 
canonical dimension. For this purpose we expand around some constant field po 

r*M = f d d x\u k (po) + U' k (po)(p - Po) + lu' k \po)(p - Po) 2 + • • • 



7 See |5^| for early applications of the derivative expansion to critical phenomena. For a recent study 
on convergence properties of the derivative expansion see |53| . 

8 See also S, M p3] for the importance of expanding around (f> — (f>o instead of <f> = and refs. |5q, pTl p8[. 
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- l - (z k {po) + Z' k (p )(p - p ) + \zl{ PQ ){p - p ) 2 + • • • ) WVa 

+l(z k (po) + z' k ( Po )(p -(*>) + . ..)r{d^y<t> a 

~Y k (p )pd^p +...}. (2.33) 

The field po may depend on k. In particular, for a potential U k with minimum at 
Po(k) > the location of the minimum can be used as one of the couplings. In this 
case po(fc) replaces the coupling U' k (p ) since U' k (p (k)) = 0. In three dimensions one 
may start by considering an approximation that takes into account only couplings with 
positive canonical mass dimension, i.e. U' k (0) with mass dimension M 2 and U k (0) with 
dimension M 1 in the symmetric regime (potential minimum at p = 0). Equivalently, 
in the spontaneously broken regime (potential minimum for p ^ 0) we may take po{k) 
and U k (po). The first correction includes then the dimensionless parameters U k (po) 
and Z k (p ). The second correction includes U k 4 \po), Z' k (p ) and Y k (p ) with mass 
dimension M _1 and so on. Already the inclusion of the dimensionless couplings gives a 
very satisfactory description of critical phenomena in three-dimensional scalar theories 



2.4 Flow equation for the average potential 

For a discussion of the ground state, its preserved or sponataneously broken symmetries and 
the mass spectrum of excitations the most important quantity is the average potential £4(p). 
In the absence of external sources the minimum p of U k ^o determines the expectation value 
of the order parameter. The symmetric phase with unbroken O(N) symmetry is realized if 
po(k —>■ 0) = 0, whereas spontaneous symmetry breaking occurs for po(k — > 0) > 0. Except 
for the wave function renormalization to be discussed later the squared particle masses M 2 
are given by M 2 ~ U'(po = 0) for the symmetric phase. Here primes denote derivatives with 
respect to p. For p ^ one finds a radial mode with M 2 ~ U'(p ) + 2p U"(p ) and N — 1 
Goldstone modes with M 2 ~ U'(po). For vanishing external sources the Goldstone modes are 
massless. 

Therefore, we want to concentrate on the flow of U k (p). The exact flow equation is obtained 
by evaluating eq. ( |2.17| ) for a constant value of <p a , say a (x) = 4>5 a i, P = \<t> 2 - O ne finds the 
exact equation || [| 

_ If d d q dR k (N-l 1 \ 

mM = 2 J wr^r {-Mo" + m ) ' (2 - 34) 

with 

M (p,q 2 ) = Z k (p,q 2 )q 2 + R k (q) + U' k (p) 

M^q 2 ) = Z k (p,q 2 )q 2 + R k (q) + U' k (p)+2pU / k \p) (2.35) 
parametrizing the (a, a) and (1,1) element of T k 2 ^ + R k (a ^ 1). 
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As expected, this equation is not closed since we need information about the p and q 2 - 
dependent wave function renormalizations Z k and Z k for the Goldstone and radial modes, 
respectively. The lowest order in the derivative expansion would take Z k = Z k = Z k (p ,k 2 ) 
independent of p and q 2 , so that only the anomalous dimension 

7 1 = -^\nZ k (2.36) 

is needed in addition to the partial differential equation ( p.34j) . For a first discussion let us 
also neglect the contribution ~ dtZ k in dtR k and write 



^-U k (p) = 2v d k d 



{N-i)in u M]^ d J u ' k{p)+2pmp) 



z k k* y u v z k k 2 



(2.37) 



with 

„-^ 2 ^r(|),„ 2 ^,„3^,, 4 ^. (2.38) 
Here we have introduced the dimensionless threshold function 

4 J (ZTT) a q 2 + Z k R k {q) + k 2 w 

It depends on the renormalized particle mass w = M 2 / (Z k k 2 ) and has the important property 
that it decays rapidly for w ^> 1. This describes the decoupling of modes with renormalized 
squared mass M 2 /Z k larger than k 2 . In consequence, only modes with mass smaller than k 
contribute to the flow. The flow equations ensure automatically the emergence of effective 
theories for the low-mass modes! The explicit form of the threshold functions ( [2.39| ) depends 
on the choice of R k . For a given explicit form of the threshold functions eq. (|2.37| ) turns into a 



nonlinear partial differential equation for a function U depending on the two variables k and 
p. This can be solved numerically by appropriate algorithms |6(J . 



Eq. ( |2.37|) was first derived || as a renormalization group improved perturbative expression 
and its intuitive form close to perturbation theory makes it very suitable for practical 
investigations. Here it is important to note that the use of the average action allows for 
the inclusion of propagator corrections (wave function renormalization effects) in a direct and 



systematic way. Extensions to more complicated scalar models or models with fermions |30 
are straightforward. In the limit of a sharp cutoff and for vanishing anomalous dimension eq. 
( [2.37| ) coincides with the Wegner-Houghton equation || for the potential, first discussed in 
(see also ©0,0). 

Eq. fl2.37| ) can be used as a practical starting point for various systematic expansions. 
For example, it is the lowest order in the derivative expansion. The next order includes 
g 2 -independent wave function renormalizations Z k (p), Z k (p) in eq. ( |2.35| ). For N = 1 the 
first-order in the derivative expansion leads therefore to coupled partial nonlinear differential 
equations for two functions U k (p) and Z k (p) depending on two variables k and p. We have 
solved these differential equations numerically and the result for different values of k is plotted 
in fig. p]. The initial values of the integration correspond to the phase with spontaneous 
symmetry breaking. More details can be found in sect. [3]. 
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Figure 1: Average potential Uk(<f>) for different scales k — e*. The shape of is displayed 
at smaller intervals At = —0.02 after the minimum has settled at a constant value. This 
demonstrates the approach to convexity in the "inner region", while the "outer region" becomes 
k-independent for k — > 0. 

3 0(iV)-symmetric scalar models 

In this section we study the iV-component scalar model with 0(iV)-symmetry in three and two 
dimensions. Due to "triviality" the four-dimensional potential flows rapidly to the perturbative 



domain and needs not to be discussed in detail here. The O(N) model serves as a prototype 
for investigations concerning the restoration of a spontaneously broken symmetry at high 
temperature. For N = 4 the model describes the scalar sector of the electroweak standard 
model in the limit of vanishing gauge and Yukawa couplings. It is also used as an effective 
model for the chiral phase transition in QCD in the limit of two quark flavors. In condensed 
matter physics N = 3 corresponds to the Heisenberg model used to describe the ferromagnetic 
phase transition. There are other applications like the helium superfluid transition (N = 2), 
liquid- vapor transition (N = 1) or statistical properties of long polymer chains (N = 0). 

3.1 Scaling form of the exact flow equation for the potential 



Let us first derive an explicitly scale-invariant form of the exact flow equation ( |2.34| ). This 
will be a useful starting point for the discussion of critical phenomena. It is convenient to use 
a dimensionless cutoff function 
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and write the flow equations as 



<9 /N f°° , d ( x \ ( N - 1 1 \ , . 



with 



Sk © = m rk © - vn (10 = - 2x h k (10 ~ Wk {&) ■ (3 - 3) 



We parametrize the wave function renormalization by 

Z fc (p,0) Z fc (p, 0) - Z fe O, 0) d _ 2 
* fc (p) = — - , z k (p ) = 1 , py k (p) = k , (3.4) 

so that 



-^ = 2^[(iV-l)/ ^^;r ?) ^j+/ d ( * ;ry,^ + Z fe py fc fc 2 - d )]+ A^, (3.5) 



2\ 



l i{w\7],z) = \ I dyyh k (y)[(z + r k (y))y + w] 1 . (3.6) 



where /q(w; r/, 2) is a generalized dimensionless threshold function (y = x/k 

poo 

2 7o 

The correction A( k contributes only in next (second) order in a derivative expansion and will 
be omitted. Finally, we may remove the explicit dependence on Z k and k by using scaling 
variables 

u k = U k k~ d , p = Z k k 2 ~ d p, K = Z k k 2 - d p (3.7) 

Evaluating the t-derivative at fixed p and using the notation u' = du/dp etc. one obtains the 
scaling form of the exact evolution equation for the average potential 

d t u Vp = -du + (d-2 + n)pu' + 2v d { (N - l)l d (u; 77, z) + l d (u' + 2pu", 77, z + py) } (3.8) 

All explicit dependence on the scale k or the wave function renormalization Z k has disappeared. 
Reparametrization invariance under field scaling is obvious in this form. For p — > a 2 p one also 
has Z k — > or 2 Z k so that p is invariant. This property needs the factor Z k in R k . This version 
is therefore most appropriate for a discussion of critical behavior. The universal features of 
the critical behavior for second-order phase transitions are related to the existence of a scaling 
solution^]. This scaling solution solves the differential equation for a /^-independent function 
u(p), which results from ( |3.8| ) by setting dtu = 0, and similar for z(p) etc.. In lowest order in 
the derivative expansion with a constant wave function renormalization the scaling potential 
can be directly obtained by solving the second-order differential equation dtu = 0. It has 
been shown that, of all possible solutions, the physical fixed point corresponds to the solution 
u(p) which is non-singular in the field [58| p^] . For d = 3 the only nontrivial solution 



corresponds to the Wilson-Fisher |7C| fixed point 



9 More precisely, d t u' — is a sufficient condition for the existence of a scaling solution. The universal 
aspects of first-order transitions are also connected to exact or approximate scaling solutions. 
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In first nontrivial order in the derivative expansion one has to supplement two partial 
differential equations || || for the p-dependent wave function renormalizations z and 

z = z + py 

d t z = r/z + pzid — 2 + 77) 

-{4v d /d)p~ l {mi fi {w,z,ri) - 2m( jl (w,w,z,z,r}) + m d t2 (w,z,r))} 
-2v d (z - z)p- 1 {l{{w, z, 77) - (2/d){z - z)l d+2 (w, ~z, 77)} 
-2v d z'{(N - l)lf(w, z, rj) - {8/d)ni A {w, w, z, z, V ) 

+ (5 + 2z"p/z') l d (w, z, rj) - (4/d)z'pl d+ j 2 (w, w, z, z, T])}, (3.9) 



d t z = rjz + pz'id — 2 + i]) 

-2v d (z' + 2pz")l((w, 5, rj) + 8v d pz'{3u" + 2pu'")l d (w, z, rj) 

+Av d (2 + l/d) p(z') 2 l d 2 +2 (w, z, rj) - (8/d)v d p(3u" + 2pu'"fm d A {w, z, rj) 

-(16/d)v d pz'(3u" + 2pu"')m d+2 (w, z, rj) - (8/d)v d p{z') 2 rh d+ \w, z, rj) 

+ (N - l)v d { -2(z'- p~\z - z)) lf(w, z, rj) - (8/d)p(u") 2 m d (w, z, rj) 

- (16 /d)pu"z'm d+2 {w, z, r]) - {8/d)p(z') 2 m d+ \w, z, rj) + 4(5 - z)u"l d {w, z, rj) 

+4 [z\z -z) + {l/d)p{z') 2 ) z, r,)} (3.10) 

They invoke the "threshold functions" t d lU2 = {l d 1>ri2 , m d 1;ri2 , m* lif , a , n d in2 } defined by the 
integrals 

CO 

c ^ = -U d ™*~ % {(p(y)+w)*m ' (3 ' n) 



with X — 1, y(d y p) 2 , y(d y p) 2 , yd y p for I, m, rh, n respectively. We have defined 
w = u', w = u' + 2pu", z = z + py, p(y) = y(z + r(y)), and p(y) = y{z + r\yj). The derivative 
d t only acts on the fc-dependence of the cutoff Rk, i.e. d t p(y) = —y(v r (y) + ^yd y r(y)) and we 
note that d t d y p = d y d t p. Finally we abbreviate l d Q = l d etc., and l d Q3.6Q corresponds to the 
rule (p + w)~ n — > — log(p + w). The expression for the anomalous dimension rj = —dtlnZk 
can be obtained from the identity dtZk(nk) = 0, which corresponds to a definition of Zk by 
<2fc(ft) = 1- (For N = 1 one has y = 0, z = z). 

These equations are valid in arbitrary dimension. We will show that they lead to 
quantitatively accurate results. This is done by a numerical solution with initial values 
specified at a microscopic scale k = A. 

3.2 Threshold functions 

In situations where the momentum dependence of the propagator can be approximated by 
a standard form of the kinetic term and is weak for the other 1PI correlation functions, the 
"non-perturbative" effects beyond one loop arise to a large extent from the threshold functions. 
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We will therefore discuss their most important properties and introduce the notation 



n + ^nfi -1 u2n-d 
4 V d 



k 2n ~ d J -0^- d d t R k (q)(Z k zq 2 + R k {q) + Z k wk 2 y^ +1 \ 



l d +1 (w;r),z) 

n + o„,o ow 

l d (w; 77) = l d (w; rj, 1) , l d (w) = l d n (w; 0, 1) , l d n = l d (0)- 







(3.12) 



The precise form of the threshold functions depends on the choice of the cutoff function R k (q). 
There are, however, a few general features which are independent of the particular scheme: 



1. For n = d/2 one has the universal property 

ll n = 1. 



(3.13) 



This guarantees the universality of the perturbative /3-functions for the quartic coupling 
in d = 4 or for the coupling in the nonlinear cx-model in d — 2. 

2. If the momentum integrals are dominated by q 2 ~ k 2 and R k (q) ^ k 2 , the threshold 
functions obey for large w 

l d {w) ~ w- {n+1) . (3.14) 

3. The threshold functions diverge for some negative value of w. This is related to the fact 
that the average potential must become convex for k — > 0. 

It is instructive to evaluate the threshold functions l d (w) explicitly for a simple cutoff 
function of the form 

R k = Z k (k 2 - q 2 )Q(k 2 - q 2 ) (3.15) 

where (x = q 2 ) 



5(x- k 2 )(k 2 -x) 



l d {w) = k 2 ~ d dxx^w + l^ + k 2 -' dxx^- , . 

J J x + (k 2 — x)(d(k 2 — x) + k 2 w 

(3.16) 

This corresponds to an "optimal cutoff" in the sense of \T7\. The second term in the 
expression for l d has to be properly defined and we consider eq. ( |3.15| ) as the limit 7 — > 00 of 
a family of cutoff functions, e.g. 



Rk 



27 
1 + 7" 



■Z k k 



6 



7 2 + l 7 2 + 2' 



q M*5 



2 \ 7 



cxp 



2 7 (q 



1 + 7 V^ 2 



2 \ 7 



(3.17) 

This confirms that the second term in eq. ( |3.16| ) vanishes and one obtains the threshold 
functions 



l d (w) 



d(l + w) 



2n 

d(l + w) n+1 



for n > 1 



(3.18) 
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In leading order in the derivative expansion (z — 1, y — 0) and neglecting the anomalous 
dimension in Iq, this choice of the threshold function yields the simple evolution equation 

3,u = -*, + (i - 2 + „)p u< + ^ + 1+M , 1 +2 ^, } • (3.19) 

In particular, for A/" = 1, d = 3 and neglecting the anomalous dimension we end with a simple 
partial differential equation 

ft»=-3 M + Pu' + ^, + (3.20) 



For our numerical calculations we mainly use the smooth exponential cutoff (|2.7|) which is 
qualitatively similar to ( |3.15|) and corresponds to (|3.17|) with 7=1. 



3.3 Second-order phase transitions 

At a second-order phase transition there is no mass scale present in the theory. In particular, 
one expects a scaling behavior of the rescaled effective average potential u k (p). Let us follow a 
typical trajectory which describes the scale dependence of the derivative u' k (p) as k is lowered 
from A to zero. Near the phase transition the trajectory spends most of the "time" t in 
the vicinity of the fc-independent scaling solution given by dtu'^(p) = 0. Only at the end 
of the running the "near-critical" trajectories deviate from the scaling solution. For k — > 
they either end up in the symmetric phase with k = and positive constant mass term m 2 
so that u' k (0) ~ m 2 /k 2 ; or they lead to a non-vanishing constant po indicating spontaneous 
symmetry breaking with k — > Z k 2 ~ d po. The equation of state involves the potential Uq(p) for 
temperatures away from the critical temperature. Its computation requires the solution for 
the running away from the critical trajectory. 

More precisely, we start at short distances (k = A) with a quartic potential or 

«a(p) = A a (p-« a ). (3.21) 

We arbitrarily choose Aa = 0.1 and fine tune k\ so that a scaling solution is approached at 
the later stages of the evolution. There is a critical value K cr for which the evolution leads to 
the scaling solution. For the results in fig. 1 a value slightly larger than K cr is used. As 
k is lowered (and t turns negative), u' k (p) deviates from its initial linear shape. Subsequently 
it evolves towards a form which is independent of k and corresponds to the scaling solution 
dtu'*{p) = 0. It spends a long "time" t - which can be rendered arbitrarily long through 
appropriate fine tuning of k\ - in the vicinity of the scaling solution. During this "time" , the 
minimum of the potential u' k (p) takes a fixed value k*, while the minimum of U k (p) evolves 
towards zero according to 

|Oo(A0 = knjZ k . (3.22) 

The longer u' k {p) stays near the scaling solution, the smaller the resulting value of po(k) when 
the system deviates from it. As this value determines the mass scale for the renormalized 
theory at k = 0, the scaling solution governs the behavior of the system very close to the phase 
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transition, where the characteristic mass scale goes to zero. Another important property of the 
"near-critical" trajectories, which spend a long "time" t near the scaling solution, is that they 
become insensitive to the details of the short- distance interactions that determine the initial 
conditions for the evolution. In particular, after u' k (p) has evolved away from its scaling form 
w*(p), its shape is independent of the choice of Aa for the classical theory. This property gives 
rise to the universal behavior near second-order phase transitions. For the solution depicted 
in fig. 1, k finally grows in such a way that po(k) approaches a constant value for k — > 0. 

As we have already mentioned the details of the renormalized theory in the vicinity of 
the phase transition are independent of the classical coupling Aa. Also the initial form of 
the potential does not have to be of the quartic form of eq. ( |3.21| ) as long as the symmetries 
are respected. The critical theory can be parameterized in terms of critical exponents |7D| . 
An example is the anomalous dimension 77, which determines the behavior of the two-point 
function at the critical temperature according to eq. (EOT) for k — > 0. Its value is given 



by the value of rj for the scaling solution. The critical exponents are universal quantities 
that depend only on the dimensionality of the system and its internal symmetries. For our 
three-dimensional theory they depend only on the value of N and can be easily extracted from 
our results. We concentrate here on the exponent v which parameterizes the behavior of the 
renormalized mass in the critical region and consider the symmetric phase 

m 2 = L^M = fcV fc (0) for k - 0. (3.23) 
Z k dp 

The behavior of m 2 in the critical region depends only on the distance from the phase 
transition, which can be expressed in terms of the difference of ka from the critical value 
n cr for which the renormalized theory has exactly m 2 = 0. The exponent v is determined from 
the relation 

m 2 ~\5k k \ 2u = \k k - k ct \ 2u . (3.24) 

Assuming proportionality 8k& ~ T c — T, this yields the critical temperature dependence of 
the correlation length £ = m~ l . For a determination of v from our results we calculate m 2 
for various values of near n cr . We subsequently plot ln(m 2 ) as a function of In \Sk^\. This 
curve becomes linear for — > and we obtain v from the constant slope. 
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N 


V 


V 





0.589' 
0.590 s 


0.5882(ll) a 
0.5875(25) b 
0.5878(6) cl 
0.5877(6) e 


0.040' 
0.039 s 


0.0284(25) a 
0.0300(50) b 


1 


0.643' 
0.6307 s 


0.6304(13) a 
0.6290(25) 6 
0.63002(23) c2 
0.6294(9) e 


0.044' 
0.0467 s 


0.0335(25) a 
0.0360(50) 6 
0.0364(4) c2 
0.0374(14) e 


2 


0.697 / 
0.666 s 


0.6703(15) a 
0.6680(35) 6 
0.67166(55) c2 
0.6721(13) e 


0.042' 
0.049 s 


0.0354(25) a 
0.0380(50) 6 
0.0381(3) c2 
0.042(2) e 


3 


0.747' 
0.704 s 


0.7073(35) a 
0.7045(55) b 
0.716(2) cl 
0.7128(14) e 


0.038' 
0.049 s 


0.0355(25) a 
0.0375(45) 6 

0.041(2) e 


4 


0.787' 
0.739 s 


0.741(6) a 
0.737(8) 6 
0.759(3) cl 
0.7525(10) e 


0.034' 
0.047 s 


0.0350(45) a 
0.036(4) 6 

0.038(l) e 


10 


0.904' 
0.881 s 


0.894(4) cl 
0.877 d 


0.019' 
0.028 s 


0.025 d 


100 


0.990' 
0.990 s 


0.989 d 


0.002' 
0.003 s 


0.003 d 



Table 1 Critical exponents v and r] for various values of N. For comparison we list results 
obtained with other methods as summarized in ]7j|/ , |7l}/ , ]7^/ , fiZ^J , fi73iJ: 

a) From perturbation series at fixed dimension including seven-loop contributions. 

b) From the e- expansion at order e 5 . 

c) From high temperature expansions (cl: /]75]/ ; c2: /|7^|/ ; see also f(7h\ 

d) From the 1/N -expansion at order 1/N 2 . 

e) From lattice Monte Carlo simulations ]7£} , [7S|/ (see also fi8(\ , \8l\]). 

f) Average action in lowest order in the derivative expansion j\Ib\ \^ J. 

g) From first order in the derivative expansion for the average action with field dependent wave 
function renormalizations (see j[L$ for N = 1 and for N > 1). 

In table 1 we compare our values for the critical exponents v and 77 obtained from 
the numerical solution of the partial differential equation ( [3.8D with results obtained from 
other methods (such as the e-expansion, perturbation series at fixed dimension, lattice high 
temperature expansions, Monte Carlo simulations and the 1/iV-expansion). We show both the 
lowest order of the derivative expansion (f) - which needs an additional equation for rj - 
and the first order (g) which corresponds to solutions of the system of eqs. (|3.8|) , ( |3.9|) , ( ft.lOQ . 
As expected rj is rather poorly determined since it is the quantity most seriously affected 



77/). 
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by the omission of the higher derivative terms in the average action. The exponent v is in 
agreement with the known results at the 1 % level, with a discrepancy between lowest and first 
order ((/) vs. (g)) roughly equal to the value of 77 for various N. Similar results are found for 



a variety of forms for the infrared cutoff function g(| [5^, ^3], [84| . We observe a convincing 
apparent convergence of the derivative expansion towards the high precision values obtained 
from the other methods. 

3.4 Kosterlitz-Thouless transition and essential scaling 



In two dimensions (d = 2) and for N = 2 the Kosterlitz-Thouless phase transition ]55[ may 
describe the critical behaviour of various two dimensional systems. It poses a challenge to 
our theoretical understanding due to several uncommon features. The low temperature phase 
exhibits a massless Goldstone boson like excitation despite the fact that the global U(l) 
symmetry is not spontaneously broken by a standard order parameter In this phase the 
critical exponents depend on the temperature. In the high temperature phase the approach 
to the transition is not governed by critical exponents but rather by essential scaling. These 
features can be explained by a description of the phase transition as a condensation of vortices 



| 85fl . Within our renormalization group equation there is no direct need for the introduction of 
vortex degrees of freedom since all nonperturbative effects are included in the exact equation. 
We therefore aim here at a complementary description in terms of the linear O(N) scalar model 
without ever introducing vortices. This constitutes an excellent test for the nonperturbative 
content of our equation. 

Already a quartic polynomial truncation for the potential shows the characteristic 



properties of the low temperature phase |p7|. This section is based on 82] and uses the 



first order in the derivative expansion ( |3.8j ), (|3.9|) , (|3~T0| ). One finds, indeed, essential scaling 



in the high temperature phase. Our discussion of the low temperature phase is expected to 
hold at the 10-20 percent accuracy level. 

In this section we present a unified picture for O(N) models in d — 2. The most important 
parameter is the location k of the potential minimum and we investigate its dependence on k. 
For k 3> 1 the evolution is dominated by the N — 1 Goldstone modes (N > 1). More precisely, 
the threshold functions at the minimum vanish rapidly for w = 2hu"(k) ^> 1. For d = 2 and 
N > 2 the effective coupling of the nonlinear a- model for the Goldstone bosons is given by 
k~ 1 ||. The universality of the (3-f unction for the nonlinear coupling |88] implies for k ^> 1 
an asymptotic form (N > 2, d = 2) 

TV 2 N — 2 

dtK = P K = — + -ttt^k- 1 + 0(k- 2 ). (3.25) 

47T l07T z 

In the linear description one can easily obtain an equation for k by using du' k (nk)/dt = 
d t u' k {n) + u'l{K)dtK = together with equation (|3.8| ). By evaluating the above equations for 
large k it is possible to compare with (|3.25|) . In the much simpler quartic truncation one 
already obtains || the correct lowest order in the above expansion. In order to reproduce the 
exact two loop result one has, however, to go even beyond the first order in the derivative 
expansion [Q . 
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From equation (|3.25| ) we expect that k will run only marginally at large k. As a consequence 
the flow of the action follows a single trajectory for large —t and can be characterized by a 
single scale. Notice that the perturbative /^-function vanishes for N = 2 since the Goldstone 
bosons are not interacting in the abelian case. Thus for large k one expects a line of fixed 
points which can be parameterized by k. This fact plays a major role in the discussion of the 
Kosterlitz-Thouless transition below. It is responsible for the temperature dependence of the 
critical exponents. We have evaluated f3 K = pjp(N - 2)/(4tt) + 0k-\N - 2)/(16tt 2 ) + ... 
numerically from the solution of eqs. (|3.8| ) - ( p.lOj ) and extracted the expansion coefficients 
for large k (see table 2). Our numerical result is close to the two-loop expression. 



N 


Cerge Cjfg C s 




3 
9 

100 


2.81±0.30 2.94 1.00 
1.22±0.03 1.25 1.05 
1.08±0.04 1.02 1.06 


1.00 0.79 
1.00 0.84 
1.00 0.87 



Table 2 Nonabelian nonlinear sigma model in d = 2. We show the ratio between the 
renormalized mass tur and the nonperturbative scale Aerge in comparison with the known 
ratio ^ invoking A^: C ERGE = m R /A ERGE , Cj^ = itir/Aj^, C s = m R /k s . We also 
display the expansion coefficients for the beta function. 

For the nonabelian nonlinear sigma model in d = 2, N > 2 there exists an exact expression 
f89fl for the ratio of the renormalized mass rriR and the scale Aj^ which characterizes the 
two loop running coupling in the MS scheme by dimensional transmutation. The flow 
equation (2.17) together with a choice of the cutoff and the initial conditions also 



defines a renormalization scheme. The corresponding parameter Aerge specifies the two 
loop perturbative value of the running coupling k _1 similar to Ajjg in the MS scheme. The 
numerical solution of the flow equation permits us to compute wir/ A ERGE . (Two loop accuracy 
would be needed for a quantitative determination of A E rg E /Aj^.) In table 2 we compare our 
results with the exact value of m R /A-j^. The quantitative agreement between C E rq E and 
Cjis is striking and suggests A ERGE rs Ajjg. We also report the ratio m R /k s with k s defined 
by n{k s ) = 0. 

The abelian case, iV = 2, is known to exhibit a special kind of phase transition which is 
usually described in terms of vortices ||85|| . The characteristics of this transition are a massive 
high temperature phase and a low temperature phase with divergent correlation length but 
zero magnetization. The anomalous dimension rj depends on T below T c and is zero above. It 
takes the exact value 1]^ = 0.25 at the transition. The most distinguishing feature is essential 
scaling for the temperature dependence of tur just above T c , 

m R ~e~(^, C=~- (3-26) 

We have already mentioned the existence of a line of fixed points for large values of k, 
which is relevant for the low temperature phase. The contribution of a massless (Goldstone) 
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cQ_ y 0.15 - 




K 



Figure 2: The beta function for d = 2, N = 2. 



boson in the renormalization group equation (w(k) = 0) is responsible for the finite value of 
r\. This in turn drives the expectation value of the unrenormalized field to zero (even for a 
nonvanishing renormalized expectation value k), 



p = K/Z k ~ K{k/\f. 



(3.27) 



We can observe this line of fixed points to a good approximation (cf. fig. 2), although the 
vanishing of f3 R is not exact (we find f3 K ~ —4 ■ 10 « + 0(k~ 2 )). Along the line of fixed 
points the anomalous dimension rj differs from zero. This reconciles the massless Goldstone 
boson with Coleman's no to theorem |9(J for a free propagator. The line of fixed points ends 
at a phase transition which corresponds to a microscopic parameter ka* as an initial value of 
the flow. In order to verify essential scaling we have to examine the flow for values of k\ just 
below that point, ka = ^a* + <5ka, Ska ~ — (T — T c ). Then crosses zero at the scale k s and 
we find the mass by continuing the flow in the symmetric regime (minimum at k — 0). In fig. 
3 we plot (In (m^/A))~ 2 against ka and find excellent agreement with the straight line ( p.2€ ). 

How does f3 K have to look like in order to yield essential scaling? Since there is only one 
independent scale near the transition, one expects mn(T) = C s k s (T), where k s denotes the 
scale at which k vanishes, i.e. n(k s ,T) = 0. For k close to and below ft* we parameterize f3 K 
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Figure 3: Essential scaling for d = 2, N 
function of k a = ^a* — H(T — T c ). 



2. The renormalized mass tur is plotted 



(this approximation is not valid for very small k 



i 

V 



K 



(3.28) 



For conventional scaling one expects ( = and the correlation length exponent is given by v. 
Integrating equation ( |3.28| ) yields for ( ^ 0, 5k = k — k*: 



ln(fc/A) 



(3.29) 



For k = k s the first term ~ (—5k)~^ is small and independent of T (since —5n(k s ) = k*) 
and equation ( ^.29| ) yields the essential scaling relation fl3.26p for ( = 1/2. Usually, the 
microscopic theory is such that one does not start immediately in the vicinity of the critical 
point and the approximation (|3.28| ) is not valid for k « A. However, if one is near the critical 
temperature the trajectories will stay close to the critical one, K c (t), with k c (0) = ka*. This 
critical trajectory converges rapidly to its asymptotic value and (3 K gets close to eq. ( |3.28 



at some scale A' < A. As a result, one may use equation (|3.29|) only in its range of validity 
(k < A') and observe that 5k\/ is also proportional to T c — T. The numerical verification of 
( P.28[ ) is quite satisfactory: Fitting our data yields = 0.248, ( = 0.502 and u^ 1 = 2.54. The 
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Figure 4: Details of the beta function for d = 2, N = 2. The curve is a fit to eq. (|3.28|) 



uncertainty for ( is approximately ±0.05. The numerical values of f3 K and the approximation 
are shown in fig. 4. 



N 


V 


V 



1 
2 


0.70 0.75 
0.92 1 


0.222 0.2083... 
0.295 0.25 
0.287 0.25 



Table 3 Critical exponents v and rj for d = 2. We compare each value with the exact result. 

One can use the information from figure [5] or [| in order to determine and therefore 
77* = f](n*), the anomalous dimension at the transition. We plot rj against n in fig. 5. One 
reads off 77* = 0.287 ± 0.007 where the error reflects the two methods used to compute and 
does not include the truncation error. For ka > ka* or T < T c the running of k essentially 
stops after a short "initial running" towards the line of fixed points. One can infer from 
figure H] the temperature dependence of the critical exponent rj for the low temperature phase. 
In summary all the relevant characteristic features of the Kosterlitz-Thouless transition are 
visible within our approach. 
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Figure 5: Temperature dependence of the anomalous dimension rj for the low temperature 
phase, d = 2, N = 2. The line of fixed points is characterized by k and ends in the critical 
point for the Kosterlitz-Thouless phase transition. We also show the flow towards the line of 
fixed points and the flow in the high temperature phase away from the critical point (left). 
The spacing between the points indicates the speed of the flow. 

We also have computed the critical exponents for the "standard" second-order phase 
transitions for N = 1 (Ising model) and N = 0. They are displayed in table 3. The 
exponents are less accurate than the ones for d = 3 displayed in table 1 (g). This confirms 
that the convergence of the derivative expansion is controlled by the size of rj. Nevertheless, 
the agreement with the exact results is better than 10 %. 

We conclude that the first order in the derivative expansion of the exact flow equation for 
the effective average action gives a quantitatively accurate picture of all phase transitions of 
scalar models in the O(N) universality class for arbitrary dimension 2 < d < 4. As the next 
step, a reliable error estimate would be very welcome. 
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